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It is suggested that marked features of symmetry breaking mechanism and elementary excitations 
in chiral helimagnet come up as visible effects in electron spin resonance (ESR) profile. Under 
the magnetic field applied parallel and perpendicular to the helical axis, elementary excitations are 
respectively described by the helimagnon associated with rotational symmetry breaking and the 
magnetic kink crystal phonon associated with translational symmetry breaking. We demonstrate 
how the ESR spectra distinguish these excitations. 



In magnetism, chirality means the left- or right- 
handedness associated with the helical order of mag- 
netic moments. Helimagnetic order can arise from spon- 
taneous symmetry breaking in systems with competing 
exchange interactions T| ( "symmetric" helimagnets), or it 
can be stabilized by the Dzyaloshinkii-Moriya (DM) an- 
tisymmetric exchange interaction [2, 0], which is real- 
ized in crystals lacking rotoinversion symmetry ("chi- 
ral" helimagnets). To clarify physical outcome of the chi- 
ral spin modulation is of great interest, especially in con- 
nection with the symmetry breaking mechanism and the 
spectrum of elementary excitations which are quite sen- 
sitive to the direction of the applied magnetic field. 

In a chiral helimagnet, under the magnetic field par- 
allel to the helical axis, the ground state (GS) generally 
changes from planar spiral to conical states [Fig. [BJa)]. 
The incommensurate modulation period 27r/(5o is fixed 
through Qo — ta,n~^{D/J), where D and J are nearest- 
neighbor DM interaction and ferromagnetic exchange in- 
teraction strengths 0, S]- The GS has infinite degen- 
eracy associated with arbitrary choice of the origin of 
the phase angle ipo. Consequently, the rotational sym- 
metry around the helical axis is spontaneously broken. 
Then, there appears helimagnetic spin-wave (chiral heli- 
magnon) mode[6] as the Nambu-Goldstone (NG) mode, 
which is well described in conventional spin wave picture. 
The chiral helimagnon has been studied in the context of 
cubic magnet MnSi[3, Q and its peculiar nature has at- 
tracted much attention in its own right 

On the other hand, under the magnetic field applied 
perpendicular to the helical axis, the GS possesses a peri- 
odic array of the commensurate (C) and incommensurate 
(IC) domains partitioned by discommensurations (DCs), 
i.e. the internal lattice which is called magnetic kink 
crystaHMKC) or sometimes referred to as chiral soliton 
latticed, Q is stabilized as shown in Fig. [T](b). Actually, 
formation of the MKC state is reported in CuB90a[1C|. 
This state is also regarded as non-trivial topological GS. 
The topological GS in chiral magnet has attracted active 
attention from various viewpoints pj|. As the magnetic 
field strength increases, the spatial period of MKC lat- 



tice, Lkink, increases and finally goes to infinity at the 
critical field strength. Recently, we showed that this in- 
ternal lattice exhibits mutual sliding which may be ex- 
perimentally detectable [l^, In this case, the GS has 
infinite degeneracy associated with arbitrary choice of 
the center of mass position. Consequently, the transla- 
tional symmetry along the helical axis is spontaneously 
broken. Then, the elementary excitations are described 
by "phonon" mode of correlated kinks. What is interest- 
ing is that we can control the size of the first Brillouin 
zone of the MKC lattice upon changing the magnetic field 
strength. 




FIG. 1: (a) Conical and (b) magnetic kink crystal (MKC) 
states. The helical axis is z-axis. 

The elementary excitations in the kink crystal state 
was first investigated by Sutherland[lJ|. He considered 
the sine-Gordon model for a single scalar field corre- 
sponding to the tangential tp-mode of the planar XY spins 
and found that the elementary excitations consist of the 
acoustic and optical bands separated by the energy gap. 
The acoustic band is formed out of correlated transla- 
tions of the individual kinks and corresponds to gapless 
NG bosons. The optical band corresponds to renormal- 
ized Klein-Gordon bosons. In chiral helimagnet, we need 
to take account of not only the (p-mode but the longitu- 
dinal 6'-mode {9 is an angle between the spin vector and 
the helical axis). In previous works we pointed out 
that the 9-mode acquires an energy gap originating from 
the DM interaction. 

Then, natural question arises as to whether the heli- 
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magnon and MKC phonon have observable consequences 
for the magnetic response using ESR technique. In this 
paper, we demonstrate how the symmetry breaking pat- 
terns and the elementary excitations come up in the ESR 
signals. 

In the ESR experiment, the static magnetic field Hq 
is applied to cause Larmor precession of magnetic spins. 
Then supplying electromagnetic energy carried by mi- 
crowave radiation, resonant absorption occurs at the pre- 
cession frequency. The microwave is described as the uni- 
form oscillating magnetic field (rf field) h(i) polarized 
in the direction perpendicular to Hq (Faraday configu- 
ration). The rf field gives rise to the Zeeman coupling 
with spin, Tiz = — H(t) • So, where H(t) — cje^J-B^it) 
[ge is the electron's g-factor and /is is the Bohr mag- 
neton) and So is the uniform {q = 0) component of 
the spin variable. For H(t) = iJ^e^ cos (tjt), the ESR 
spectrum (absorbed energy per unit time) is given by, 
Q{oj) — ujH'j^x'li.iii^) /2, where e'' with /x = x,y,z de- 
notes the unit vector along x, y, and z axis [Fig. 1], re- 
spectively, and w is a microwave frequency. The imag- 
inary part of the dynamical susceptibility, x'/ii^ {^) — 
(l - e.-^^^^^) {uj) /2, is related to the correlation 
function C^v (w) = {Sq [uj) Sq) through the fluctuation- 
dissipation theorem. In quantum mechanical language, 
the Lamor precession corresponds to equally spaced Zee- 
man splitting of the energy levels. Because of the equal 
spacing of the quantum energy levels, the quantum- 
classical correspondence exactly holds and the classical 
frequency is equal to quantum one as far as we consider 
Gaussian fluctuations. 

First, we consider the case where the magnetic field 
is applied parallel to the helical axis (z-axis) and the 
rf field is polarized along the y-axis. Then, the ele- 
mentary excitations are described as spin waves over 
the conical magnetic structure. A quantized spin wave 
is helimagnon. Then, the ESR spectrum is given by 
Qhmag (w) = ujH^Xyy {^) /2- To compute Xyy (w), we as- 
sume that the magnetic atoms form a three dimensional 
lattice and a uniform ferromagnetic coupling exists be- 
tween the adjacent chains to stabilize the long-range or- 
der. Then, the Hamiltonian is interpreted as an effective 
one-dimensional model based on the interchain mean field 
picture and is written as, 

j 

- j5]5|5|+i + /f^5](5|)2 - Ho • ^S„ (1) 

j j j 

where Sj represents a spin located at the j-th site along 
the helical axis (z-axis) and — S'^ ± iS^. The mono- 
axial DM vector is D = De^ and J = \J + iD\ = 
\/\P + D^. The lattice constant is c. We include the 
easy-plane anisotropy with strength K±. For Hq = 0, 




FIG. 2: (a) Helimagnon dispersions for Ho/Hoc ~ 0, 0.7, and 
1.. We took D/J = 0.5 and K±/J = 2. Black dots indicate 
the location of the resonance energies, (b) Field dependence 
of the resonance energy as a function of Hq/Hqc- 

the planar helical structure is stable under the condition 
K±^ / J > 1 — y^l + [D / JY which is assumed to be sat- 
isfied. For < Ho < Hoc = 2S{J - J + K^), the GS 
is described by Sf = 5'e±*('3o^3+¥'o) sin 6', where the cone 

angle is given hy 6 = 6o ^ cos-^[Hol {2S{j - J + K^)}]. 

To obtain the spin wave spectrum, we rotate the basis 
frame of the crystal coordinate {e+,e~,e^} to the basis 
frame of the local coordinate {e^,e~,e|} where the di- 
rection of e| points to the equilibrium spin direction at 
the j-th site[15|. In the spirit of conventional spin- wave 
approximation, we obtain the spectrum, 

^ = ^[l-cos(9c)] [A-7C0s(gc)], (2) 

where g is a wave number, A = 1 + {K±/J) sin^ 6o and 
7 — (J/J) sin^ 9o + cos^ 9o- This result reduces to the 
one obtained by Kataoka 7] and Maleyev[8] using the 
continuum approximation (q limit). In Fig. ^a), 
we show the helimagnon dispersion for Ho = 0, 0.7Hoc, 
and Hoc- Upon increasing the field, linear dispersions for 
< Ho < Hoc continuously crosses over to the quadratic 
dispersion hujq = 2JS'(1 — cosg) at Ho — Hoc- The 
Goldstone mode at g = corresponds to the rigid ro- 
tation of the whole helix. For Ho > Hoc, the equilib- 
rium state is forced-ferromagnetic state and the spin wave 
spectrum acquire the field-induced gap. 

Now, it is straightforward to obtain the helimagnon 
resonance spectrum, 

ttS 

Qhmag (uj) = -^Wff,^(5 [uj - WqJ 
^ [(4o + "Qo)' + ^oKo - "Qo)']' (3) 

where u^^ = \/{P/^Qo ± 1) /2 and P = S\2J + 
K±_ sin^ 00 - J{1 + cos^ Oq + (J/J) sin^ Bo}]. Note that 
the external uniform field couples to the q = ±Qo com- 
ponent of the spin wave excitation, since the field is seen 
in the local frame as spatially rotating field with modu- 
lation wave- number Qo- Consequently, we have a single 
branch of resonance energy, as shown in Fig.[2I^b). As we 
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shall see, this situation drastically changes in the case of 
the MKC phonon resonance. 

ESR signal in chiral helimagnet MnSi was reported by 
Date et aL16]. At that time, however, they adopted the 
formula obtained by Yoshimori^l|| and Cooper et al.fv^ 
for symmetric helimagnetic structure stabilized by frus- 
tration among the exchange interactions^] . In the case of 
symmetric helimagnet, the spin wave dispersion exhibits 
dips at g = ±Qo and the corresponding energy gaps van- 
ish for K± = 0(13. There are no such additional dips 
in chiral helimagnon spectrum. We see, however, it may 
not be easy to distinguish the spin wave spectra of chi- 
ral helimagnet from those of symmetric helimagnet sim- 
ply by ESR profile, because both cases give apparently 
quite similar field dependence of the resonance energies 
as shown in Fig. EJb). 

Next, we consider the MKC phonon resonance when 
the magnetic field is applied perpendicular to the heli- 
cal axis (y-axis) and the rf field is polarized along the 
z-axis. The MKC state is described in terms of the 
slowly varying polar angles 6{z) and f{z). The vector 
spin density is defined by S (z) — J2j ~ ^j) — 

(sin cos sin6'(z) sin (/3(z), cos6'(z)) . Then, mini- 
mizing the continuum version of the Hamiltonian ([T]) , we 
obtain the MKC state as a stationary state described 
by 9 = tt/2 and cos[(Pq{z)/2] ~ sn(2iirz/Lkink), where 
Lkink — 8KE/ttQo is the period of the MKC lattice. K 
and E denote the elliptic integrals of the first and second 
kind, respectively, with the elliptic modulus k {0 < n < 
1). "sn"is Jacobi-sn function. The IC to C transition 
occurs at the critical field strength H^/JS = {ttQq/A) 
at which Lkink diverges. The elliptic modulus k is deter- 
mined by the condition ^ Hq /Hq — k/E (k) . The IC to 
C transition in chiral magnet is actually reported in real 
materials [ill . For example, in the case of Cri /^NbS9[l9t. 
iJg takes values from about 1 to 1.4 kOe and in the case 
of CUB2O4 20], from 0.5 to lOkOc depending on temper- 
atures. 

In this case, the rf field couples with S^{z,t) = 
S cos 6{z,t) and ESR spectrum is given by Qph (t^) = 
uoHlx'L (^) /2- To compute x'L {'^)^ need the explicit 
form of the propagating mode S'^(z, t) ~ —Su{z, t) where 
u{z, t) = 6{z, t) — 7r/2 describes small fiuctuation around 
the MKC state. Although full description should include 
the (^-mode, the rf field couples to only 0-mode and it is 
enough to consider the 0-mode only. By using the mode 
expansion for u{z,t), we set up the vibrational Hamilto- 
nian given as collections of harmonic oscillators fl5]. The 
explicit form of the quantized phonon wave function is 
given by 



q n— — oo 



,^-i(q-nGmKc)z+iujgt-^'\ ^ ^ 



(4) 

where (hq) are the phonon creation (annihilation) op- 
erators. The crystal-momentum q and the eigenfre- 



quency ojq are expressed in terms of a real parameter 
a running over —K' < a < K' , where K' is the complete 
elliptic integral of the first kind with the complementary 
modulus k' . For the acoustic branch, q = lQ^[Z{a, k!) -t- 
-Ka/iiKK')] and hujq = £o^A2-f k'2 sn2(a,K')/2. For 
the optical branch, q = lQ^[Z{a^ k') -\- ira/ (2KK') + 
dn (a, k')cs (a, k')] and hujq — £o\/A2 -|- sn"^ (a, n') /2. 
Z(a, k') is the elliptic zeta-function. We here intro- 
duced the characteristic length and energy units /q — 
Lkink/2if = 4£;/7rQo - Qo ^ and Eq = JS^c/Iq = 
JS^ttQqc/AE ~ DS^c, respectively. It is essential that 
the energy gap A = y^SE/n — 2 opens at g = 0, because 
of the existence of the DM interaction. The Fourier coef- 
ficients Un can be computed by performing contour inte- 
gral of the real space wave function given in Ref . [13] . To 
obtain the ESR spectrum, we need ^ it/ {2K^/k!) and 
Un = -i-^d'{^^di{iTTan/2K)/ smh['Kan/{2K)-nTiK' /K] 
for n ^ [an is determined by the resonance condition 
^ given below], 'di is the Jacobi theta function and 
^1 = ^i(O)- Since di{iTTa/2K) is purely imaginary, all 
Un are real. 

The reciprocal lattice constant of the MKC lattice is 
given by 
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Lkink ^KE 



(5) 



The first Brillouin zone of the MKC lattice is ]g] < 
Gmkc/2 and the energy gap between the acoustic and 
optical branches opens at the zone boundary. As limiting 
forms, we have hiOq ~ et^^j^^^^rv^ (v is constant) for 
q < Gmkc/2 and hug ~ £0 \qlo\/V2 for Gmkc/2 < q. 
Now, we are ready to understand the ESR by the MKC 
phonon. Since the rf field along the z-axis carries the 
wave number g = 0, the resonant absorption is caused 
by the MKC phonon modes with a series of special wave 
numbers 



q = qn = ?lG; 



MKC- 



(6) 



The correlation function can be easily computed by using 
Eq. Q and we obtain the ESR absorption spectrum. 



n=0 



-(5 (w - cj„) , 



(7) 



where lOn — ■ This expression together with C/„ and 
bJn complete a closed formula for the MKC phonon res- 
onance. For n = 0, the bottom of the acoustic branch 
(q = and a = 0) gives = For n > 1, the optical 
branch contributes to the resonance. 

As the magnetic field increases from zero to , Gmkc 
decreases from Qo to zero. On the other hand, the origi- 
nal atomic lattice constant c gives natural cutoff and the 
atomic Brillouin zone boundary ±27r/c irrespective of the 
external magnetic field. Usually, 211 /Qq ~ 10c— 100c and 
therefore Gmkc is much smaller than the atomic zone 



4 




(b) 







Q 


Q 


(c-1) 

day 








J 


I 








1.0 

/ 


{ ' 


2.0 

hay 

"e7 



asymptotic form of the resonance frequencies for large n, 




0.6 O.S 1 

Hn/Hf, 



(c-2) 

day 



So 



FIG. 3: Energy dispersion of the MKC phonon in the reduced 
zone scheme for (a-1) smaller and (a-2) larger magnetic field 
strengths. The vertical broken lines indicate the Brillouin 
zone boundaries q — ±Gmkc/2. (b) Resonance energy ujn for 
n = to n = 10 as functions of Hq/Hq. We took D/J = 0.5 
and K^l J = 0. The derivative absorption dQp\^/duj for (c-1) 
Ho/Ho = 0.8 and (c-2) Ho/Hq = l-10-^ In (a-1), (a-2), (c- 
1), and (c-2), black dots indicate the location of the resonance 
energies. 



boundary 27r/c. In Fig.[3Ua), we schematically depict that 
the distribution of the resonance energy levels becomes 
more and more dense upon increasing the magnetic field 
strength. In Fig. Elb), we show the resonance energies 
LOn for 71 = to n = 10 as functions of Hq/Hq. To obtain 
the result presented here, we first numerically solve Eq. 
(pl in terms of the parameter a, and then compute the 
corresponding frequency w„, which completes evaluation 
ofEq. (O. In Fig. [3J c-1), we show the derivative absorp- 
tion dQph (w) /duj for Hq/Hq = 0.8. The delta function 
is replaced by 6{lli) = 7r^^e/{Lj^ + e^) with e = 10^'*. Al- 
though the weight \Un\^ rapidly decays for higher order 
resonances, the peak structure becomes visible by taking 
the derivative. 

Of special interest is the region in the vicinity of the 
IC-C transition, where the distribution of the resonance 
levels are quite dense. In Fig. [3lJc-2), we show the 
case for Hq/Hq 1. We see that a series of many 
densely spaced resonance lines appear. Using the rela- 
tion K{k) ~ log(4/vT^-«?) and E{k) ~ 1 which hold 
for K < 1, we have k ~ y/ Hq /Hq and therefore obtain the 



£o 



K 



log (4/^1 - Hq/H* 



(8) 



Conversely, a series of resonance fields for large n are 
given by HQn/HQ ~ 1 — 16 exp (— 27rneo/Sa') for a fixed 
frequency ijJ. Our energy unit Eq ~ DS^c usually 
amounts to J/100 — J/10, corresponding to ImeV to 
lOmeV in energy scales. These energy scales correspond 
to microwave frequencies in THz region (quantitative de- 
tail depends on Eq). So, our effects should be detectable 
in submillimeter wave ESR measurements jllj. 

We stress that the MKC phonon resonance never 
occurs in the symmetric helimagnet due to energetic 
frustration|lj, where not the MKC but the "fan" struc- 
ture is stabilized under the field perpendicular to the he- 
lical axis Physical background behind this dif- 
ference is that in chiral helimagnet the crystallographic 
chirality plays a role of "topological protectorate" for the 
MKC lattice state to appear as the stable GS. 

In basic physical ideas, the effects we proposed here is 
one of examples to detect spin dynamics of phase mod- 
ulated states by polarized probes such as inelastic neu- 
tron beam or X-ray. For example, neutron beams can 
probe the MKC phonon mode via the differential cross 
section (Pa/dVlduj oc (1 — k^) (S'^(cj)S'ij^) where lj and k 
are respectively frequency and scattering wave-number of 
the neutron[24|. Then, the scattering event occurs when 
both momentum conservation, = q — nGiviKCi ^^'^ 
the energy conservation, lu = iwq, are satisfied, where q 
is the MKC phonon wave-number. The polarized X-ray 
beam may also detect the MKC state via the generalized 
spin-orbit coupling between the spin magnetic moment 
and X-ray. These topics will be treated separately in a 
forthcoming paper. 

Finally, we make theoretical comments on the MKC 
state. The MKC apparently seems to be one-dimensional 
object which is fragile against three dimensional cou- 
plings. However, it is not necessary to worry about this. 
Many features in physics of incommensurate magnets 
may be understood based on the Ginzburg-Landau free 
energy with a non-uniform order parameter as a function 
of f/iree-dimensional coordinates. Then one should select 
a solution minimizing the functional that corresponds to 
the modulated phase. As a result of this analysis, we 
find a structure with a modulation along one axis in the 
crystal is easily stabilized. In such a case, it is enough to 
take into account the invariant involving derivatives with 
respect to one coordinate (z-coordinate in the present 
case). More rigorously speaking, we need to exclude a 
possibility that we have a structure with multiple mod- 
ulation vectors in a single crystallographic domain. But 
it is known that the realization of this kind of structure 
is hard to occur (see for example, Ref.^]). This is the 
reason why we can safely start with the effective one- 
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dimensional model as we did in this paper. 
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